Analysis of the resonance peak and magnetic coherence 
seen in inelastic neutron scattering of cuprate superconductors: 
a consistent picture with tunneling and conductivity data 
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Assuming the exchange of antiferromagnetic spin fluctuations as the Cooper pairing mechanism 
we calculate the doping dependence of the resonance peak seen in inelastic neutron scattering and 
the magnetic coherence effect. Most importantly, we find that the resonance peak in the magnetic 
susceptibility, Im x(ci,<^), appears only in the superconducting state, that it scales with Tc, and that 
magnetic coherence is a result of a d-wave order parameter. We further analyze the structure of 
Im X below Tc, the position of the peak at Ures and the consequences for photoemission, tunneling 
spectroscopy and the optical conductivity. 



74.20. Mn, 74.72.-h, 74.25.-q 

For analyzing the pairing mechanism in high- Tc su- 
perconductors it is important to understand the spin- 
excitation spectrum as observed by inelastic neutron 
scattering (INS) In particular the doping and tem- 

perature dependence of the spin susceptibihty Imx(q, uj) 
and its relationship to the superconducting transition 
temperature Tc are important. INS experiments show 
the appearance of a resonance peak at w^es only be- 
low Tc [|l| and find a constant ratio of uJres/Tc — 5.4 
for underdoped YBa2Cu307_5 (YBCO) and overdoped 
BiaSrzCaCuaOg+i (BSCCO) Recent INS data on 

La2_a;Sr3;Cu04 (LSCO) reveal strong momentum- and 
frequency-dependent changes of Imx(q, a;) in the super- 
conducting state which the authors called magnetic co- 
herence effect Q. In particular, Imx(Qi) with Qj — 
(1 ± (5, 1 ± 5)ii is strongly suppressed compared to its 
normal state value below uj meV, while it increases 
above this frequency. Furthermore, the incommensurate 
peaks become sharper in the superconducting state 
Several theoretical approaches have been considered for 
the resonance peak its relation to the spectral 



function seen in ARPES and tunneling data [|12 13|, and 
for the magnetic coherence fl^ , but no unified theory 
has been presented. Here, using an electronic theory, we 
perform calculations to demonstrate the significance of 
the feedback of superconductivity on Im x- This permits 
us to derive for underdoped and overdoped cuprates the 
relationship between different measurements like optical 
conductivity and tunneling. So far the resonance peak 
and magnetic coherence have been treated separately and 
in both cases without taking into account this important 
feedback of superconductivity on Im x- 

In this Communication we use an electronic theory for 
the spin susceptibility and for the Cooper pairing via 
exchange of antiferromagnetic spin fluctuations to ana- 
lyze the consequences of the superconducting feedback on 
magnetic coherence and the resonance peak, and on the 
relationship between INS, tunneling, and optical conduc- 



tivity. The investigation of the effect of superconductiv- 
ity on Im X is an extension of our previous work |]l0| , |l5| . 
Using RPA and self-consistent FLEX |16| calculations 
for Imx(q, w), we present results for the kinematic gap 
(or spin gap) o^o, ujres, ^res/Tc, and the gap function 
A(a;) in reasonable agreement with experiments. Most 
importantly, we find that our electronic theory can ex- 
plain consistently within the same picture inelastic neu- 
tron scattering, optical conductivity, and SIN tunneling 
data. Furthermore, the same physical picture gives re- 
sults for underdoped and overdoped cuprates . 

In order to analyze the kinematic gap and the position 
of the resonance peak it is instructive to start with the 
bare Lindhard BCS susceptibility Q at q = Q = (tt, tt) 



Imxo(Q, c^) = i 5] {(1 - 2f{Ek)) d{uj + 2Ek) 



+ {2,f{Ek)-l)5{uj-2Ek)} 



(1) 



where f{Ek) is the Fermi function and E\^ = -^/e^ + is 
the dispersion of the Cooper-pairs in the superconducting 
state. We use a gap function with d-wave symmetry, 
Ak = Aq (cos kx — cos ky)/2, which can be calculated self- 
consistently within our FLEX-approach. For the normal 
state dispersion, we employ a tight-binding band 

ek = — 2t [cos kx + cos ky — 2t' cos k^ cos fcy — /i/2] . (2) 

Here, t is the nearest neighbor hopping energy, t' denotes 
the ratio of next-nearest neighbor to nearest neighbor 
hopping energy, and /i is the chemical potential. We use 
t' as a fitting parameter in order to describe the Fermi 
surface topology of both materials YBCO and LSCO. 
In Eq. (|^) we take t' — Q and, for simplicity, we do 
not consider a bilayer coupling via a hopping integral t j_ 
Jl9| , p0t . Imxo(Q,t^) involves two characteristic frequen- 
cies. The first, ludoSi arises from the density of states 
of the Bogoluibov quasiparticles (i.e. the Cooper-pairs) 
which have a gap in their spectrum due to superconduc- 
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FIG. 1. Numerical results for the resonance peak and mag- 
netic coiierence in the weak-coupling limit, (a) Imaginary part 
of the RPA spin susceptibility (in units of states/eV) versus 
u) in the superconducting state at wavevector q = Q = (tt, tt) 
for U/t = 1, 2, 3, and 4 (from bottom to top). As in Ref. 
[0 we find ujre-a ~ 41 meV. Below the kinematic gap loq, 
Im x(Q,'^) is zero, (b) Calculated magnetic coherence: the 
solid curves correspond to the superconducting state whereas 
the dotted curve is calculated in the normal state. The ob- 
served four peaks occur at Qi = (1 ± 5, 1 ± 5)-k and in the 
figure we show only the peaks at Q; = (1, 1 ± 5)-n. In our cal- 
culations we find 5 ~ 0.18. These results are in fair agreement 
with experiments, see Ref. [^]. 

tivity, ludos — 2A(x,T). Here, x is the doping concen- 
tration. The second, coq, at which Imxo(Q:'^) starts to 
increase represents the existence of a d-wave supercon- 
ducting order parameter and is the so-called kinematic 
gap [[7y2Cf|. We point out that, for i' > 0.3 the kinematic 
gap is washed out. 

In Fig. |l| (a) results for the spin susceptibility 



Imx(Q,t^) 



(1 - C/Rexo(Q,c^))2 + C/2 (Imxo(Q,^)= 



(3) 



are shown . Here, U stands for an effective Hubbard in- 
teraction. We choose t' = 0.2 in order to describe the 
Fermi surface topology of both YBCO and LSCO. Gen- 
erally, one finds that the structure of Im x is determined 
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FIG. 2. Consistent picture of INS and tunneling data, 
(a) Imaginary part of the RPA spin susceptibility at 
q = Q = (7r,7r) calculated within the FLEX approxima- 
tion. For the normal state (dotted) we get Usf — O.lt and 
in the superconducting state we obtain ujrea = 0.15t. As- 
suming t = 250 meV we find that 0.16t — 40 meV. Inset: 
Imaginary part of the gap function at T = 0.075r for wave 
vector q ~ (7r,0). (b) Calculated density of states for the 
same parameters and temperatures as in (a). 



by Im xo if {URe xo) 1 and by {URe xo) = 1 if this 
can be fulfilled. We again find the two characteristic fre- 
quencies luq and LUres — ^DOS which Imx is peaked. 
Furthermore, one clearly sees that for increasing U the 
peak in Im x shifts to lower energies, and most impor- 
tantly, becomes resonant for U = Ucr which satisfies the 
condition 



1 



(4) 



which signals the occurence of a spin-density-wave col- 
lective mode. The real part is given (at T = 0) by: 



E]^E]^+Q, — CkCk-l-Q — AkAk-i-Q i5k 



k 



2i!^k£'k+Q 



(5) 



Re Xo(Q:i^res) has been investigated in detail in Ref. 
]lO{ , where it was found that the spin-density- wave col- 
lective mode, that satisfies Eq. (||), can explain the dip 
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FIG. 3. Calculated FLEX results ujres versus doping. In 
the overdoped regime, with Tc oc Aq we find a constant 
ratio of uires/Tc — 8. 



and hump feature observed in the photoemission spec- 
tra on BSCCO |^l|. In particular, it was shown that 
the broad humps are at the same position for both the 
normal and superconducting state. 

We find from Eq. (|^) that, in the normal state 
where no resonance appears, the spin wave spectrum is 
mainly determined by the spin fluctuation frequency uj^f 
(roughly the peak position) and 



Imx(Q,t^ = 2 I 2 



(6) 



In the superconducting state one finds that Im x peaks 
resonantly at uj^es where tOres — 2A as it can be already 
seen from Eq. (^. More precisely, we find for opti- 
mal doping, where Eq. (Q) determines the structure of 
Im X the important relation ujres{T) = 2Ao(T) —UsfiT). 
Physically speaking, the resonance peak appearing in INS 
only below Tc is mainly determined by the maximum of 
the superconducting gap, but renormalized by normal 
state spin excitations. This provides a simple explana- 
tion for the observed 41 meV resonance peak in optimally 
doped YBCO |l[], because Raman data suggest 2 A ~ 58 
meV 



and LOsf ~ 17 meV (at 100 K) as extracted from 
NMR experiments 

In Fig. |l|(b) we show results for the q-dependence of 
Imx(q, i^)- We perform our calculations for U — 2t and 
a superconducting gap of 2 A = 10 meV as measured in 
Raman scattering in optimally doped Lai.85Sro.i5Cu04 
p^ . For u = 10 meV we get two peaks at q = Qi. In 
the superconducting state we find a sharpening of the 
peaks due to the occurrence of a gap. This simply means 
that the lifetime of the quasiparticles is enhanced due 
to a reduced scattering rate. At 4 meV these peaks are 
strongly suppressed as seen in the experiment Q. Fur- 
thermore, we find no signal for u; < 4 meV. This is due 
to the kinematic gap seen in Fig. ^(a) which is indepen- 
dent of q. Note, the situation would be totally different 
if LSCO would have an isotropic gap where all states for 




FIG. 4. Calculated scattering rate in the normal state for 
T = 1.5Tc (dotted curve) using equation (^) and supercon- 
ducting state 0.75Tc (solid curve) using Kubo formula and 
yielding a threshold (2A-|-ci;res). The results are in fair agree- 
ment with Ref. Ei. 



< w < 2Ao ~ 20 meV would be forbidden. In this case 
no kinematic gap (or spin gap) would be observed. Thus, 
we can conclude from Fig. m that, already in the weak 
coupling limit where no lifetime of the Cooper-pairs (i.e., 
A independent of lu) is considered, we are able to explain 
the resonance peak and the magnetic coherence effect. 

In order to consider the important feedback effect 
of A on the spin excitation spectrum, we now dis- 
cuss our results obtained in the strong-coupling limit 
(i.e. w-dependent) solving self-consistently the general- 
ized Eliashberg equations within the FLEX approxima- 
tion Jl5| , p^ . Note that, only U/t and the tight-binding 
dispersion relation e(k) (with its band filling fi) enter the 
theory as free parameters. 

In Fig. ^(a) we present results for Imx(Q,w) calcu- 
lated for U = At and for an optimum doping concentra- 
tion X = 0.15 {fi = 1.65). In the normal state (dotted 
curve) we find cusj ~ O.lt, whereas for T < Tc the res- 
onance peak (solid curve) appears at LOres = 0.15i. The 
dashed curve corresponds to T = 0.9Tc where the super- 
conducting gap starts to open. Thus, the peak position 
reveals information on the temperature dependence of 
the superconducting gap. For temperatures T < 0.75Tc 
the resonance peak remains at LUrcs — 0.1 5i and only the 
peak height increases further. We find that the height 
of the peak is of the order of the quasiparticle lifetime 
l/T{LUres), whcrc r(k, a;) = Lului Z(k, a;)/Re Z(k, lj). Z 
denotes the mass renormalization within the Eliashberg 
theory. Thus we can conclude that the resonance peak 
becomes observable because the scattering rate decreases 
drastically below Tc [|lO|. 

In Fig. ||(b) we show the corresponding calculated den- 
sity of states N{lj). Below T < 0.75rc we find that the 
value of 2A determined from peak-to-peak stays approx- 
imately constant and is very close to the value ujres seen 
in INS, i.e. 41 meV, as shown in (a). This is in good 



3 



agreement with measured SIN tunneling data in Ref. [g5j . 
However, in SIN tunneling a renormalized value of 2A is 
observed. Note, a direct measurement of A(w) (e.g. SIS 
tunneling) would lead to higher values. For example, we 
show in the inset of Fig. |^ (a) the imaginary part of gap 
function at wave vector q ~ (tt, 0) where the gap has its 
maximum. It is peaked at w = 0.25t. 

In Fig. 1^ we show results for ujres as a function of 
the doping concentration. We find that, for a fixed U 
Eq. (|^) cannot be fulfilled in the overdoped case [ p6[ . 
Thus we find that in this regime the resonance peak is 
mainly determined by Im xo(Qj and thus by 2Ao. On 
general grounds one expects Tc oc Aq in the overdoped 
regime, where the system behaves mean-field (BCS) like. 
Recently, this has been confirmed within the FLEX ap- 
proach [|l8|. Thus, we conclude that ujres/Tc should be 
a constant ratio. We find ujres/Tc — 8 which is larger 
than the observed value in BSCCO ||] . This is due to an 
underestimation of Tc within FLEX. 

In contrast to the overdoped case we find in the under- 
doped regime, where Tc oc (ris denotes the superfluid 
density) , that the resonance condition Eq. (||) yields 
ujres LOsf which is decreasing. Note that, the super- 
conducting gap guarantees that Eq.(||) is fulfilled. Thus 
we find a decreasing resonance frequency for decreasing 
doping in agreement with earlier calculations |ll|,0 . To 
summarize our discussion we have the following result: 




underdoped 
jgf, optimal doping 
overdoped 



(7) 



where the optimally doped case corresponds to Xopt — 
0.15 holes per copper site. This predicted doping depen- 
dence of the resonance peak position should be further 
tested experimentally. 

In order to discuss the consequences of our analysis in 
particular the feedback of superconductivity on Im x for 
various superconducting properties we derive ||2^] 



[/2 

Im Efoj) = 

4 



dio' 



coth — - 
2T 



tanh 



2T 



(8) 



where N{llj) — X^k'^d^l ~ -^k) is the density of states. 
Eq. is valid in both the normal and superconducting 
state. It permits discussion of how much S(u;) reflects 
ujef and LVres, for example. We see that the feedback of 
superconductivity on Im x causes approximately a shift 
for the elementary excitations lu —^ lu + Aq for the super- 
conducting state in the spectral density Im x{Q,^) /4. 
Using Eq. in Eq. (|^) would not take into account the 
important feedback of superconductivity on Im x- Eq- 
can be used to demonstrate the relationship between 
INS and optical conductivity measurements for the nor- 
mal state, but note for the superconducting state we cal- 
culate a from (GG -I- FF). Using Drude theory we find 



that the scattering rate t~^{uj) agrees qualitatively with 
— 2Im S(w) for the normal state. However, in order to 
get a quantitative agreement with experimental data one 
has to use t~^{uj) = r(Q,aj). This is shown in Fig. ||, 
in a good agreement with Ref. From this analysis 

we can conclude that optical conductivity data and thus 
the lifetime of the quasiparticles, as well as the resonance 
peak and SIN tunneling data, can be understood within 
our electronic theory. 

In summary, we are able to explain consistently all 
characteristic facts about the spin excitation spectrum 
and its doping dependence in high- Tc cuprates within an 
electronic theory. In particular we find that the reso- 
nance peak is a rearrangement of spectral weight of the 
normal state which happens only below Tc- Furthermore 
we show that magnetic coherence is connected with the 
resonance peak and can be explained by a kinematic gap 
and by a d-wave symmetry of the superconducting or- 
der parameter. By taking into account the feedback of 
superconductivity on Imx(q, w) we argue that ARPES 
results, tunneling data, and measurements of the opti- 
cal conductivity are consistent. This further strengthens 
spin fluctuation exchange as the relevant mechanism for 
Cooper-pairing in high- Tc cuprates. 
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